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Abstract. We present a consistent derivation of the recently proposed 5D anisotropic 
standing wave braneworld generated by gravity coupled to a phantom-like scalar field. 
We explicitly solve the corresponding junction conditions, a fact that enables us to 
give a physical interpretation to the anisotropic energy-momentum tensor components 
on the brane. So matter on the brane represents an oscillating fluid which emits 
anisotropic waves into the bulk. We also analyze the Sturm-Liouville problem 
associated to the correct localization condition of the transverse to the brane metric 
and scalar fields. It is shown that this condition restricts the physically meaningful 
space of solutions for the localization of the fluctuations of the model. 
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1. Introduction 

Since the early proposals of braneworld models involving large extra dimensions and 
4D delta-function sources with both positive and negative tensions [H |2] there has been 
a lot of activity in this area with the aim of solving several open questions in modern 
physics (see [31 H] for reviews). Most of these models were realized as time independent 
field configurations. However, quite soon, mostly within cosmological approaches, there 
have appeared several braneworlds that assumed time-dependent metrics and fields 
in an attempt to address several open problems in astrophysics and cosmology [5]- 
[12]. In general, these braneworld models have a series of remarkable features: they 
develop a novel geometrical mechanism of dimensional reduction based on a curved 
extra dimension, provide a realization of the AdS/CFT correspondence to lowest order, 
take into account the self-gravity on the brane through its tension, can trap various 
matter fields (except for gauge bosons) inside the brane, recast the hierarchy problem 
into a higher dimensional viewpoint, and also lead to cosmological backgrounds with 
consistent dynamics that recover general relativity results under suitable restrictions on 
their parameters. All these facts provide a complex but interesting interplay between 
gravity, particle physics and geometry [1]. 

However, there are still several open relevant questions which deserve attention: 
construction of the simplest realistic solution for an astrophysical black hole on the brane 
and study its physical properties like staticity and Hawking radiation, development 
of realistic approximation schemes and numerical codes to study cosmological 
perturbations on all scales, computation of the CMB anisotropies and large scale 
structure to confront their predictions with high-precision observations 

One of the main drawbacks of the delta-function braneworld models from the 
gravitational viewpoint is related to their singularities at the positions of the branes, 
a fact that gets worse in models with more than six dimensions due to the stronger 
self-gravity of the brane. In an attempt to heal this problem, several authors smooth 
out the brane tensions [13] or introduced single [H] or several scalars jlS], and phantom 
[I6] or tachyon [17] fields as sources (for a recent review see [H]), and even considered 
modified gravity braneworlds [19]. However, scalar field thick brane configurations with 
4D Poincare symmetry also develop naked singularities at the boundaries of the bulk 
manifold if a mass gap is required to be present in the graviton spectrum of Kaluza-Klein 
fluctuations (see [HJEO], for instance). An alternative model without such a drawback 
is provided by a recently proposed de Sitter braneworld model purely generated by 4D 
and 5D cosmological constants [21]. 

In this sense, it is important to consider new generalizations that attempt to 
make more realistic the original braneworld models, or explore other aspects of higher- 
dimensional gravity which are not probed or approached by these simple models. 

In this paper we shall consider the braneworld generated by 5D anisotropic standing 
gravitational waves coupled to a phantom-like scalar field in the bulk as it was recently 
proposed in [22]. It turns out that standing wave configurations can provide a natural 
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alternative mechanism for localizing 4D gravity as well as for trapping matter fields (for 
scalar and fermion fields see [23] and [23], respectively), including gauge bosons [25] , 
which usually are not localized on thin braneworlds. A peculiarity of this anisotropic 
braneworld scenario is that it possesses non-stationary metric coefficients in the 4D part 
of the line element. We further impose Z2-symmetry along the extra dimension, a fact 
that gives rise to the need of introducing as well an anisotropic brane energy-momentum 
tensor for the self-consistency of the model. When the corresponding junction conditions 
are solved, this fact allows us to give a physical interpretation to the components of the 
anisotropic energy-momentum tensor of the brane. Then the localization of transverse 
scalar and metric fields are treated through the associated Sturm-Liouville method that 
restricts the physically meaningful parameter space of the solution. Some final remarks 
are given at the end of the paper. 



2. The model 



In this section we briefly recall the 5D standing wave braneworld model proposed in 
[22], which is generated by gravity coupled to a non self-interacting scalar phantom-like 
field [T6| [T8] which depends on time and propagates in the bulk, and is given by the 
action: 



^ (i?-2A5) + ^ W)' 



where G5 and A5 are 5D Newton and cosmological constants, respectively. It is worth 
noticing that in order to avoid the well-known problems of stability which occur with 
ghost fields, the phantom-like scalar field does not couple to ordinary matter in the 
model. 

The Einstein equations for the action ([T]) read: 

R^lu — ^SfiyR = ^T^G^T^y — k^g^y , (2) 

where Greek indices run from to 5, labeling the 5D space-time, and T^y is the energy- 
momentum tensor for the scalar field 

= -d^ct>dyct> + ]^g,ydpCt>dP(t> . (3) 

Following [22], we use the anisotropic metric ansatz: 

ds^ = e^^^'-) {dt' - e"(*'^)rfx2 - e"(*'^)rfi/2 _ ^-2u(t,r) ^^2^^ _ ^^2 ^ 

where Latin indices stand for the 4D space-time coordinates {x^ = t, ,with z = 1, 2, 3, 
are the spatial coordinates x, y, z, respectively) and the extra dimension is denoted as 
^5 _ This metric generalizes straightforwardly the thin brane metric ansatz [^, where 
A[r) ~ |r|, which is recovered in the limit when u{t,r) vanishes. The warp factor here 
A[r) is an arbitrary function of the extra coordinate r and, in principle, may model a 
thick brane configuration in the spirit of [13]. 

This braneworld constitutes a generalization of the thin brane model with the 
peculiarity that now the brane possesses anisotropic oscillations on it, which send a 
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wave into the bulk (as in [2S]), i-e. the brane is warped along the spatial coordinates 
a;, y, z through the factor e"'-*'^-*, depending on time t and the extra coordinate r. Several 
anisotropic braneworld models have been previously considered in the literature [6]-[TT] 
when addressing relevant cosmological issues like anisotropy dissipation during inflation 
[8] , braneworld isotropization with the aid of magnetic flelds [10] and localization of test 
particles [11]. Moreover, as a general feature it has been established that anisotropic 
metrics on the brane prevent the bulk from being static [9l [TO] . 

The phantom-like scalar fleld r) obeys the Klein-Gordon equation on the 
background space-time given by (j4]), 

^ d^i^/^gg^'d^cl)) = 00 = e'^^^^^c^ - 0" - 4A'0' = , (5) 

where overdots mean time derivatives and primes stand for derivatives with respect to 
the extra coordinate r. 

We further write the Einstein equations in the form: 

R^lu = -df^crdya + ^fi'^^.As , (6) 

where the gravitational constant has been absorbed in the definition of the scalar field: 

a = , (7) 

and we have used the reduced form of the energy-momentum tensor of the phantom-like 
scalar field 

T^,y = T^u - \^9t,vT = -V^aV^a , (8) 
and T = gf^-'T^,. 

It turns out that from the 05-component of the Einstein equations (|6]), it follows 
that the fields a and u are related by [22] , 

^{t,r) = ^u{t,r) , (9) 

up to an additive constant that can be absorbed into a redefinition of the 4D spatial 
coordinates in (jl]). 

On the other hand, by combining the 22- and 33-components of the Einstein 
equations ([6]), and comparing the result to the corresponding 55-component, one gets a 
quite strong restriction on the function A{r): 

A" = , (10) 

whose solution is linear in r, but corresponds to divergent warp factors. Another 
interesting fact that follows from the same analysis is that this condition also implies a 
sort of fine tuning, where the 5D cosmological constant is related to A{r) as follows: 

A5 = QA'^ ■ (11) 

It is worth noticing that the relation remains valid even if we include a self-interaction 
potential for the scalar field. 



An anisotropic standing wave braneworld and associated Sturm-Liouville problem 5 



In order to get an asymptotically vanishing warp factor and to study a more 
interesting braneworld with the above {A" = 0) behavior on the bulk, we modify the 
initial model by imposing Z2-symmetry along the extra dimension and by introducing a 
thin brane at r = 0, which, for consistency with (jl]), must be supported by an anisotropic 
energy-momentum tensor. This implies that the warp factor of the metric (jll) adopts 
the known form, similar to [2]: 

ds^ = e^"!^-! {dt^ - e'^dx^ - e'^dy^ - e-^'^dz^) - dr^ , (12) 

where a is an arbitrary constant. The non-zero components of the Ricci tensor for this 
metric read: 

3 



_-g-2a|r|^2 ^ 4^2 ^ 2a5{r) 



R-xx 


= R = e 


Rzz 


g2a r — 2u 




= - 

2 


Rrt 


3. , 

= uu 

2 



2a\r\+u 



^e-2"l^'lM - - 2ae{r)u' - 2a5{r) - \u" 



2 

[-e-^''\''\u-Aa^ + Aae{r)u'-2a5{r) + u"] , (13) 
- Aa^ - 8a5(r) , 



where e(r) is the sign function. The presence of terms proportional to 5{r) in 
the Ricci tensor, together with the anisotropy of the metric, forces us to define an 
anisotropic energy-momentum tensor on the brane r/f with different stresses along 
different directions. The simplest way to accomplish this is by proposing: 

tI: = 5{r) diag[Ao, Ai, A2, A3, 0] , Ai = A2 . (14) 

The above energy-momentum tensor can be interpreted as certain 'anisotropic effective 
fluid' which constitutes a mixture of a vacuum fluid (which is characterized just by 
a brane tension) and an anisotropic matter fluid (see |1] and references therein for a 
complete decomposition of the energy- momentum tensor induced on the brane). 

The parameter Aq is the energy density of the 'effective fluid' and Am ('^ = 1, 2, 3) 
are the stresses along the m directions (see below). In general these quantities depend 
only on time. A similar energy-momentum tensor was implemented in the context 
of Eotvos braneworlds with variable (time-dependent) tensions in analogy with fluid 
membranes in |12] . 

By taking into account the above considerations, the field equations (E]) change as 
follows: 

2 

R^lv = -d^ad^a + -g^^yk^ + SvrGsf^^ , (15) 
where the reduced energy-momentum tensor, 

V = ^A'!^ - \^9tiur , (16) 
corresponds to the matter content on the brane and takes the form: 

f^, = ^6{r) diag [(2Ao - 2Ai - A3), (Aq - Ai + \^)e\ (17) 



An anisotropic standing wave braneworld and associated Sturm-Liouville problem 6 



(Ao - Ai + A3)e", (Ao + 2Ai - 2A3)e"^", (Aq + 2Ai + A 



By making use of (fT2l) . (fT5l) and (fT7|) one can show that the relation (|9]) between the 
scalar field and the metric function is satisfied as well. The relation (fTTl) computed with 
the aid of the metric f|T2|) adopts the following form: 

As = . (18) 

Therefore, the system of Einstein and field equations reduces to a single ordinary 
differential equation, 

e-2a|r| ^ _ ^" _ 4^ae{r)u' = . (19) 

In addition to the latter, the existence of the thin brane in our setup gives rise to a 
relationship between the jump of the first derivative of the function m at r = (denoted 
by [u']) the stresses Am and the parameter a. A way to derive these relationships consists 
in integrating (fT5|) with respect to r in the range [—e, e] and then making e tend to zero, 
yielding the following junction conditions for the stresses Am^ 
4nG, 



3 



(2Ao - 2Ai - A3) 



M = (-Ao + Ai - A3) - 4a , 

[u'] = ^ ( Ao + 2Ai - 2A3) + 2a , (20) 

TlG 

a = ^ (Ao + 2Ai + A3) . 

Although the above system seems to be overdetermined since we have four equations 
for just three variables Aq, Ai and A3, it is possible to show that only three equations 
are independent. 

We then construct a standing wave solution to ( IT9l) by implementing the ansatz: 

u{t,r) = Csm{ujt)f{r) , (21) 

where C and u are some constants. In [22] it was argued that a standing wave 
configuration in the bulk would in principle provide an alternative mechanism for 
localizing gravity and binding matter fields (through an anisotropic force which drives 
the matter fields/particles toward the nodes), including light (see [23]- [25] for recent 
results). 

The equation for the radial function /(r) from (12T|) reads, 

/" + 4ae(r)/' + w2g-2aH^^Q ^ ^22) 

The general solution of this equation adopts the following form: 



/(r) = e 



-2a\r\ 



(23) 



where A, B are arbitrary constants and J2 and I2 are 2*^^ order Bessel functions of the 
first and second kind, respectively. It is worth noticing that although the Y2 Bessel 
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function is singular as r — > oo for a > 0, the whole function f l23p vanishes in this limit 
since its behaviour is dominated by the e~^"'^' factor. 

As pointed out in [22], along with the solution ( 123|) . we impose the ghost-like field 
to be unobservable on the position of the thin brane, since it oscillates together with 
the gravitational field (see ([9])). A way to accomplish this consists in setting a boundary 
condition that nullifies the ghost-like field at the position of the brane r = 0, a condition 
that quantizes the oscillation frequency, u, of the standing wave: 

^ = Xn, (24) 

\a\ 

where X„ is the n^^ zero of the 2"*^ order Bessel function J2 or Y2 depending on whether 
one takes ^ = or i3 = in ( 123|) . since the zeros of J2 and Y2 do not coincide. For an 
increasing (a > 0)/decreasing (a < 0) warp factor, the ghost-like field cr(t,r) and the 
metric function u{t,r) vanish at a finite/infinite number of points along the extra 
dimension in the bulk (nodes), forming 4D space-time 'islands' at these nodes, where 
the matter particles are assumed to be bound. 

The induced 4D Einstein space-time is locally AdS4 (in a vicinity of the positions 
of the island universes along the extra dimension) if: 

^Rab = (y{rm)qab , 0,6 = 0,1,2,3 (25) 

where "^Rab and qab are the induced 4D Ricci and metric tensors, respectively, and 
tt('^m) > is a constant depending on the position of the island r^. 

By computing the 4D induced metric on the brane islands we can write down the 
components of the 4D induced Ricci tensor as functions of the induced metric tensor and 
its first two derivatives, i.e. as functions of u{t,r) and its first two partial derivatives 
with respect to time and the extra coordinate. By looking for a proportionality between 
the induced metric and Ricci tensors on the island universes, i.e. when ti(t, r^) = 
(however, the derivatives of u at do not vanish), we see that the AdS^ effective 
character of the 4D island universes can be reached only in two distinct cases: 

i) by restricting the amplitude of the anisotropic oscillations to be very small, more 
precisely when AC << 1, for the islands which are located at finite Vm- 

ii) asymptotically along the fifth dimension (r^ — > oo). 

Thus, even when at the points where the island universes are located we have 
u(t, r^) = and the metric f|T2|) reduces to the standard 5D thin braneworld metric of [2] , 
whereas the setup ([1]) simplifies to an action describing 5D gravity with a cosmological 
constant, plus a delta-function brane as in [2], leading to the known Randall-Sundrum 
braneworld which is an AdS^ slice, in general, the metric in the islands (i.e. for each 
particular node r^) is not locally ^^5*4, but just in the above referred two cases. 

Returning to the problem of junction conditions (l20l) . if we use (|2T1) and fl23l) to 
solve them, we obtain: 

3a 



An 



5 



ACuj e{a) sm{ujt) 3a 
= 16^, ^-^-^ ~ -^'^ '-/l^l " 4^ 
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e(a) sin(a;t) , , , 3a 

'8^, - Jl) |./|a| - 



ACue{a) sin(a;t) 3a 
As = ^-7^ (-'3 - Jl) U/\a\ - -r—^ , (26) 



[m'] = — - — (Ai - A3) = ACu e(a) sm{ut) (J3 - Ji) . 

From these relations it can be seen that Aq = p is constant, and it can be interpreted 
as the energy density of the 'effective fluid' as mentioned above, while the stresses 
oscillate in the x, y and z directions with the same frequency as the function u does, 
as it was expected since it is precisely along these directions that the metric coefficients 
do oscillate in the bulk. On the other hand, the spatial components of the energy- 
momentum tensor on the brane can be consistently interpreted following a systematic 
covariant analysis developed from the viewpoint of a brane-bound observer [271 [28] . 
Let us decompose the tension components as follows: 

Ai = -p + 7r, , (2 = 1,2,3) (27) 

where the quantity p represents a constant isotropic pressure and the tTj are the 
components of the anisotropic stress that oscillate along the spatial direction x*: 
3a 



P 



47iG 



5 



ACu e{a) sinjut) 

7^1 = 7^2 = —-t; (J3 - Jl) U/|a| , (28) 

IottGs 

ACu e(a) sinjut) 

T^3= - 27ri = — {J3 - Jl) \^/la\ . 

By taking into account that in the solution for the junctions conditions (|26|) the 
parameter u}/\a\ is a zero of J2, but not of Ji and J3, in general the anisotropic stresses 
are different from zero. Furthermore, the stresses along the directions x and y are equal, 
whereas the stress in the 2;— direction tts is twice in amplitude and possesses opposite 
phase with respect to tti. This is a direct consequence of the symmetry of the metric 
under the exchange x < — > y, on the one hand, and the differences in amplitude and 
phase in the argument of the exponential metric coefficients, on the other hand. 

Let us compare the isotropic and the anisotropic parts of (fT4|) . In order to do that 
we shall consider the following physically interesting limits: 

a) The first case is Ci;/|a| ~ 1, and depending on the amplitudes A and C in the 
expressions for the components Aj, the dominant term can be the oscillatory or the 
constant one (3a/47rG5), or both terms can be of the same order. 

la) If aj/|a| ~ 1 and AC ^ 1 the constant term dominates and Aj ~ Aq = 
—Sa/AirG^. Then the oscillatory terms in Aj can be neglected compared to the 
contribution given by the constant terms and the thin brane becomes "isotropic", 
recovering the standard thin brane case [2], where the state equation is p = —p. Thus, 
the energy-momentum tensor on the brane only has a relevant contribution coming from 
the tension of the brane. 

2a) If a;/|a| ^ 1 and AC ^ 1 the oscillatory terms are dominant in Aj and 
max(Aj) ^ Ao, implying that the energy- momentum tensor of the thin brane describes 
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a kind of exotic matter given by the following expression: 



tI^ ^ 5{r) diag[Ao, Ai + Aq, A2 + Aq, A3 + Aq, 0] , Ai = A; 



'2 , 



where 



A 



ACuj e(a) sin(u;t) 



(-^3 - Ji) U/|ai, 



A3 



ACoj e(a) sin(ci;t) 



(^3 — Jl) |aj/|a| 



(29) 



SttGs 



In this case, the matter is highly anisotropic because the amplitude of the stationary 
wave is large. 

3a) If uj/\a\ ~ 1 and AC ~ 1, the oscillatory and constant terms are of the same 
order and A, ~ Aq ~ a. Due to this fact, in this case there is no way to highlight an 
interesting physical situation, in contrast to cases la) and 2a). 

b) The second interesting limit takes place when the stationary wave that propagates 
in the bulk oscillates in the high frequency regime 3> 1); in this case we have 

again three situations: 

lb) If ^ 1 and AC ^ 1 with ACXn Ci, where Ci is a finite constant, we 
have the same result as in the case 3a). Thus, the oscillatory and constant terms are 
of the same order. 

2b) If u)/\a\ ^ 1 and AC ^ 1 we have the same result as in the case 2a): the 
matter fluid is highly anisotropic. But in contrast to 2a) in this case it is because the 
frequency of the stationary wave is large, therefore, it contributes to the growth of the 
anisotropic part of the energy-momentum tensor. 

3b) If a;/|a| 3> 1 and AC ^ 1 we have the same result as in 2a): an oscillatory 
energy-momentum tensor of the thin brane which describes a kind of exotic matter. In 
this case this fact takes place because the frequency and amplitude of the stationary 
wave are large. 

Note that since u}/\a\ is a zero of J2 the quantity u}/\a\ > 1 and the third case 
u}/\a\ ^ 1 is excluded from our analysis. 

In summary, the two most relevant limits are: 

I) when the matter on the brane hardly oscillates due to the small amplitude of the 
anisotropies compared to the brane's tension, leading to a quasi-isotropic brane. This 
limit mimics the Randall-Sundrum model. 

II) when the matter on the brane is highly anisotropic and a stationary wave 
propagates through the bulk. In this case we are far from the Randall-Sundrum limit. 

It would be interesting to study a braneworld isotropization mechanism for our 
brane universe similar to those proposed in [8], HO]. In the latter model, the brane 
anisotropic energy leaks into the bulk as it evolves, a phenomenon that can also be 
interpreted from a completely 4D point of view in the framework of the AdS/CFT 
correspondence as particle production in the CFT, where energy is drawn from the 
anisotropy to fuel the process, leading to an isotropic brane over time. In principle, this 
study could be carried out for our anisotropic braneworld model since we have completely 
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solved the 5D Einstein equations, a necessary requirement for the computation of the 
projection of the Weyl tensor into the brane. 

3. Localization of transverse scalar and metric fields 

In the original paper [22] it was shown that near the nodes r^, where u{r,t) ~ 0, 
the metric f|T2|) adopts the usual thin brane form of [2] and one recovers localized 4D 
gravity on the island universes in the usual way through a massless mode solution for 
the equation that governs the dynamics of transverse traceless metric fluctuations. 

Here we shall consider the localization of small perturbations of a real massless 
scalar field defined by the action: 

S^ = -\ I Vadx^'dr g^^d.^d,^ . (30) 



2 „ 

We consider a massless scalar field since it is important to get an oscillating solution 
in our model, otherwise there will be no standing wave since the solution for the scalar 
field if will be dissipative (as in the case of including a mass term, for instance). It 
is necessary to mention that unlike [22] we study the localization of scalar field in the 
system formed by the thin brane located in r = and all the island universes. This 
will be seen below from the fact that the norm of fluctuations is defined on the whole 
domain of the extra dimension r, thus, if it is finite for the entire domain, it will also 
be finite for single islands or a set of island-universes. 

The scalar field equation corresponding to the action fl30|) is: 

^ - e-"((/^.. + cpyy) - e2>,, - e^'^l^l [4ae(r)^' + cp"] = , (31) 

where the subscripts xx, yy and zz denote second differentiation with respect to x, y 
and z, respectively. Let us propose the following ansatz for the scalar field, 

= L{r,t)P^{x)P2{,y)Pz{z) , (32) 

which transforms the previous equation into, 

'L_ _/^£l\ l£PA_e^d^_ ,,^,^ [Aae{r)L' + L"] 



Let us also set 

d'P 



(^ = 1,2,3.) (34) 

The parameters ki can be physically regarded as the momenta along the i direction of 
the brane without considering the anisotropic character of the background metric. 
By taking into account ( l34|) . we get: 

L - e^^l'-l [Aae{r)L' + L"] = - [{kl + A;2)e-" + fc^e^"] L . (35) 

In the limit u —> this equation has the well known solution: 

e'"'' , kl -k^ = 0, (36) 
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corresponding to the 4D massless scalar mode, which can be locahzed on the thin brane 
located at r = by the warp factor e^"'*^'' if a < (see, for example, |29]). 

In this paper we want to consider a simple case of transverse fluctuations to the 
island universes, i.e. when fcj ^ 0. In this case the equation (1351) takes the form: 

L - e^'^l'^l [Aae{r)L' + L"] = . (37) 
By proposing a similar to fl2T]) oscillatory ansatz: 

L{r,t) = Ksm{nt)g{r) , (38) 
where fl and K are some constants, we obtain, 

/ + 4ae(r)^' + fi2e-2'^IH^ = , (39) 

which is the same equation as ( l22l) for the metric function /(r), and thus, possesses the 
solution: 



a / V \cl\ 



(40) 



where D, E are integration constants. In contrast with the frequency of the ghost-like 
field, cu, the parameter Q in fHOj) is not quantized since in general Q ^ u. 

In order to analyze the proper way of normalization of g{r), we rewrite (139|) as: 

_ (e4aH^/y ^^2g2aM^ _ (41) 

Now we recall the Sturm-Liouville method associating to the equation, 

- [Pir)y']' + q{r)y = \s{r)y , (42) 
of a variable y{r), the norm, 

\y{r)\^s{r)dr , (43) 

where A is an eigenvalue parameter, while b and c are arbitrary real constants. By 
comparing ( 142|) with (HT]) we see that y{r) = g{r), q{r) = 0, s{r) = e^'^'''', A = fi^, 
p{r) = e^"'^' and the correct norm for (!39|) is given by. 



c 

2„2a\r 



g{r)\'e"'^''^dr , (44) 

where 6 = and c = oo. Thus g{r) belongs to the Hilbert space Hs consisting of all 
such functions for which (jH]) is finite |30] . 

Another way of looking at this problem consists of considering the action for the 
scalar field ( 130|) and taking into account ([4]) and the expression for y^, in order to show 
that its nontrivial part (the kinetic term) is, 

~ y dx^dre^^^'-^L^ + • • • , (45) 

where the dots denote the 5D contribution of the scalar field. If we further make use of 
the relation fl38|) . then the normalization condition along the extra coordinate reads: 

S^- ! dr\g{r)\\^'^\^\ , (46) 
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which precisely coincides with relation f l44p . Thus, if we want to have a localized 
transverse 5D scalar field on the brane (and due to relation ([9]), the gravitational waves 
u as well) the integral over r in ( 146|) must be finite. 

It should be pointed out here that the finiteness of the above obtained norm (l44ll 
is exactly the same for the phantom-like scalar field (0, or a), which indeed, vanishes at 
the island branes. Therefore, the latter field will be localized as well if the expression 
fj46|l constitutes a finite quantity. 

By considering the change of variable v = f2e~"'^'/|a| for the case a > we have 
the following integrating limits vi = Q/\a\ and V2 = 0. Alternatively, for a < 0, we get 
vi = fi/|a|, while V2 = 00. 

In the language of this new coordinate, the integral fH6l) adopts the form: 



where X denotes the complex conjugate of X. It is easy to get the behavior of the 
integral (1471) by direct computation. It turns out that only in the case E = and a > 
this integral is finite. Therefore the massless scalar field ip is localized on the brane at 
r = 0, as well as in all the island branes of the model since in this latter case one should 
just change the integration limits Vi and f 2, getting a finite value for each island-universe 
or for a set of them. 

It should be mentioned that if one compares our result to the result obtained in 
[22| [29] , it seems that they are in an apparent contradiction, because the scalar field is 
localized for a > in our case, while in their papers this field is localized for a < 0. 
However, the point here is that we are investigating the transverse modes and not the 
modes along the brane. Also in our work 7^ in the ansatz ( 1381) . otherwise L = 
and the scalar field would be trivial. A way of avoiding this trouble consists in 
extending the ansatz ( !38l) to contain a nonzero constant phase a 7^ 2'Kn with n ^ Z m 
the argument of the sinus function, i.e. 



In this case the = value corresponds to the lowest frequency of oscillation of the 
nontrivial mode of the scalar field which indeed adopts a constant value along the 
fifth dimension (as well as g). Then this field is localized for a < like in [22| [29] 
according to the normalization condition ( 14^ . 

Now let us roughly analyze the gravitational sector from this localizing point of 
view. The Ricci scalar for the metric ( IT2l) reads: 




+ DE J2{v)Y2{v) + \E\^ \Y2{v)\^ dv , 



(47) 



L = Ksm{Qt + a) g{r) . 



(48) 




(49) 



Thus, the gravitational 5D action is: 




(50) 
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where dots denote 5D contributions. Due to the relation ^ between the phantom-hke 
scalar field a and the gravitational field u, it turns out (as expected) that localization 
properties of these fields are similar. Thus the gravitational field u is also localized if 
we choose E = Q and a > in the solution ( 140|) under the ansatz (!38|) . 

4. Final Remarks 

In this paper we presented a consistent derivation of the 5D anisotropic standing 
wave braneworld proposed in [22] by initially assuming a Z2-symmetric factor and 
introducing a simple anisotropic energy-momentum tensor with different stresses along 
different space-time directions on the 3-brane. We also derived and explicitly solved the 
corresponding junction conditions for the braneworld model and obtained analytical 
expressions for the (oscillating) stresses of the 3-brane along the 4D spatial directions. 

By carefully looking at the energy-momentum tensor of the 3-brane proposed in 
( !T4|) we can infer from (126|) that it corresponds to an anisotropic effective fluid which is 
a mixture of a vacuum fluid characterized by its tension and an anisotropic oscillating 
matter fluid. We can also conclude that the 4D matter corresponding to this source 
has an exotic nature because it violates the weak energy condition for increasing warp 
factors and satisfles it for decreasing warp factors (a < 0). It is worth mentioning 
here that this condition is also violated in the thin brane models [2]. Besides, it 
violates as well the dominant energy condition since it is an oscillating fluid which 
emits anisotropic waves into the bulk with different amplitudes (which change sign over 
time and eventually disappear) and phases along the directions x, y and z. Nevertheless, 
one can look for an alternative mechanism that could lead to a more physical picture 
and eventually heal this drawback related to the anisotropy of the metric ansatz ( |T2l) . 
like proposing a more involved energy-momentum tensor or add some matter flelds on 
the 3-brane. An interesting situation takes place when the amplitude of the anisotropies 
is small with respect to the tension on the brane (quasi-isotropic limit), since in this 
case the braneworld is effectively isotropic. Moreover, within the framework of the 
model presented here it is possible to study the braneworld isotropization in which 
anisotropies dissipate via inflation |8] or leakage of thermal graviton radiation into the 
bulk [To], a relevant phenomenon from the brane cosmological viewpoint that deserves 
more attention. 

We showed as well the correct way of deflning the norm for the transverse scalar 
fleld ip through the Sturm-Liouville method, a fact that leads to the localization of 
this fleld just when a > and ii^ = in the solution f HOj) . However, it should be 
pointed out that this localization approach involves all the island brane universes as 
well as the thin brane located at r = 0. While the approach taken in [22] considers the 
gravity localization on each of the island universes and is related to the massless modes 
of the transverse traceless metric fluctuations of the system, i.e. to the 4D gravitons 
that live in each 3-brane. Our results seem to agree with recent results obtained in 
where localization of massive test particles about a thick braneworld with a time 
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dependent extra dimension arises just for increasing warp factors due to an oscillatory 
behaviour in time-like geodesies. Further work must be done in order to get localization 
mechanisms of matter fields with decreasing warp factors, since these are useful in solving 
the hierarchy problem. We also would like to point out that further investigations 
towards the localization of three generations of fermion fields in this model, both at the 
thin brane located at r = and at the island universes located at r^, are in progress 
and will be published elsewhere. 

We finally would like to point out that in order to avoid the instabilities related to 
the ghost-like scalar field of this braneworld we can appeal to an alternative approach 
of interpreting the phantom-like scalar field (f) as the geometrical scalar field of a 5D 
Weyl integrable manifold |311 132] , where a scalar appears through the definition of the 
covariant derivative of the metric tensor, 

D^gap = gai3d^4> ■ (51) 

This is a generalization of the Riemannian geometry, in which the covariant derivative 
of the metric tensor obeys the metricity condition, i.e. it vanishes 

D^gap = . (52) 

On the other hand, this relation indicates that the Weylian affine connections are not 
metric compatible since they also involve the scalar field in their definition: 

r^. = - I {<P,,K + - g,.r) , (53) 

where {^^} are the Christoffel symbols. As a consequence, in an integrable Weyl 
manifold specified by the pair {g^i,,(j)), the non-metricity condition (I5T]) implies that 
the length of a vector is altered by parallel transport. 
The 5D Weyl action can be written as 

= [ ^^l^e-i^ \r + 3i (V0)^ - 2^/(0)1 , (54) 

where is a Weylian integrable manifold, ^ is an arbitrary coupling parameter and 
U{(j)) is a self-interaction potential for the scalar field 0. From the formalism itself it 
becomes clear that this action is of pure geometrical nature since the scalar field that 
couples non-minimally to gravity is precisely the scalar that enters the definition of 
the affine connections of the Weyl manifold (!53|) and the non-metricity condition (15T!) 
and, thus, cannot be neglected at all in our setup. 
By performing the conformal transformation, 

diMu = e~'^giiu , (55) 
we map the Weylian action flM|) into the Einstein frame: 

r d^x^/\g\ 

where all hatted magnitudes and operators are defined in the Riemann manifold, 
^ = ^ - 1, f/(0) = e-«^f/(0) = As and we have set U{(f)) = A^e'^ in (El]) in order to 



^5^ 



(56) 
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get a 5D cosmological constant in as in ([T]). In this frame we have a theory which 
describes 5D gravity minimally coupled to a scalar field plus a cosmological constant, the 
affine connections become the Christoffel symbols, the metricity condition is recovered 
and Weyl's scalar field in ( 156|1 imitates a massless scalar field (either an ordinary scalar 
or a ghost-like scalar depending on the sign of ^ [3ll 132]). 

In the alternative approach mentioned above we can start with the action fl5^ 
which does correspond to a conventional scalar field in the Einstein frame's action fl56l) 
under the conformal transformation ( 155|) . Then we can compute the corresponding 
field equations under the metric (fT2|) (which will be different from (fTS!) since the affine 
connections involve the scalar field) and solve them in the Weyl frame. We further 
can set (p^^ = on the brane in the solution of the field equations (which also will be 
different from the solution in the Riemann manifold). Thus, in this way we can recover 
an effective 4D universe in the Einstein frame at the brane (and possibly at the island 
universes if our solution possess the same structure as the braneworld obtained in [22] ) 
completely free of ghost instabilities since the Weyl scalar corresponds to a conventional 
real scalar field. This scenario still has to be investigated in full detail, a direction which 
is currently under research. 
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